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Abstract-  This  paper  is  devoted  on  the  characterization  method 
of  RF/digital  PCB  interconnections  for  the  prediction  of  the  high¬ 
speed  signal  transient  responses.  The  introduced  method  is  based 
on  the  use  of  the  interconnection  line  RLCG-model.  Theoretical 
formulae  enabling  the  extraction  of  the  electrical  per-unit  length 
parameters  R,  L ,  C  and  G  in  function  of  the  interconnection  line 
physical  characteristics  (width,  length,  metal  conductivity,  dielectric 
permittivity  . . . )  are  established.  Then,  by  considering  the  second 
order  approximation  of  the  interconnection  RLCG-model  transfer 
matrix,  the  calculation  process  of  the  transient  responses  from  the 
interconnection  system  transfer  function  is  originally  established. 
To  demonstrate  the  relevance  of  the  proposed  model,  microwave- 
digital  interconnection  structure  comprised  of  millimetre  microstrip 
line  driven  and  loaded  by  logic  gates  which  are  respectively  modelled  by 
their  input  and  output  impedances  was  considered.  Then,  comparisons 
between  the  SPICE-computation  results  and  those  obtained  from  the 
proposed  analytical  model  implemented  in  Matlab  were  made.  As 
results,  by  considering  a  periodical  square  microwave- digital  excitation 
signal  with  2  Gbits/s  rate,  transient  responses  which  are  very  well- 
correlated  to  the  SPICE-results  and  showing  the  degradation  of  the 
tested  signal  fidelity  are  observed.  The  numerical  computations 
confirm  that  the  proposed  modelling  method  enables  also  to  evaluate 
accurately  the  transient  signal  parameters  as  the  rise- /fall-times  and 
the  50%  propagation  delay  in  very  less  computation  time.  For 
this  reason,  this  analytical-numerical  modelling  method  is  potentially 
interesting  for  the  analysis  of  the  signal  integrity  which  propagates 
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in  the  high-speed  complex  interconnection  systems  as  the  clock  tree 
distribution  networks.  In  the  continuation  of  this  work,  we  would  like 
to  apply  the  proposed  modelling  process  for  the  enhancement  of  signal 
quality  degraded  by  the  RF / digital  circuit  board  interconnection  where 
the  signal  delays  and  losses  became  considerably  critical. 

1.  INTRODUCTION 

The  last  century  was  particularly  remarkable  with  the  spectacular 
progress  of  the  microelectronic  semiconductor  industries  which  was 
till  now,  unique  in  the  mankind  history.  These  constant  technological 
progresses  are  nowadays  a  source  of  innovative  product  developments 
in  numerous  civil  sectors  as  the  mobile  phones,  multimedia  systems, 
medical  equipments  and  even  vehicles.  As  an  undeniable  quantification 
of  this  industrial  development,  according  to  Moore’s  law  prediction, 
at  each  generation  the  size  of  electronic  printed  circuit  board  (PCB) 
is  reduced  of  30%  and  every  two  or  three  years,  the  number  of  the 
transistors  integrated  in  the  microelectronic  chips  must  be  increased 
of  four  times  [1].  Despite  this  fascinating  scenario  of  technological 
development,  the  structures  of  the  current  interconnection  circuitry 
between  the  logic  gates,  electronic  chips  and  digital  devices,  and  even 
the  buses  linking  the  different  electronic  boards  as  seen  in  Fig.  1  become 
more  and  more  complicated  [2-7]. 


Figure  1.  High  speed  digital  PCB. 

Furthermore,  due  to  the  incessant  increase  of  the  operating 
frequencies  and  the  integration  density,  the  contributions  of  the 
interconnection  electric  parameters  as  the  resistance,  capacitance 
and  inductance  effects  [8-15]  cannot  presently  be  neglected  by  the 
electronic  equipment  manufacturers.  Indeed,  the  signal  integrity  which 
represents  the  nunrerical-/digital-data  is  not  preserved  especially  for 
the  high-rate  communication.  This  signal  degradation  effect  needs 
to  be  taken  into  account  by  the  electronic  designers  during  the 
manufacture  processes  of  the  high-speed  electronic  devices.  For  this 
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reason,  different  characterization  techniques  of  the  interconnection 
transmission  line  (TL)  based  on  the  5-parameter  extractions  were 
introduced  since  two  decades  ago  [16-21].  In  this  optic,  performance 
optimization  methods  dedicated  to  the  clock  tree  distribution  networks 
were  proposed  [22-24]  in  order  to  minimize  the  undesired  influences  of 
interconnection  systems  notably  for  the  VLSI  circuits.  In  addition,  a 
compensation  technique  based  on  the  insertion  of  repeaters  between 
the  piece  of  TLs  was  also  proposed  [11,25]  in  order  to  reduce 
the  TL  interconnection  negative  effects.  As  improvement  of  this 
interconnection  compensation  technique,  an  alternative  solution  based 
on  the  use  of  negative  group  delay  (NGD)  circuits  is  introduced 
in  [26,27]. 

Notwithstanding  these  various  technical  solutions,  further  simpler, 
faster  and  more  accurate  methods  need  to  be  developed  to  model 
the  interconnection  line  effects  especially  for  the  complex  circuitry 
especially,  in  transient  domain.  One  underlines  that  most  of 
existing  characterization  methods  were  developed  only  in  frequency 
domain  by  using  the  5-paranreter  characterization  deduced  from 
the  geometrical  parameters  (width,  length,  substrate  height)  of  the 
TL  and  the  electromagnetic  properties  (dielectric  permittivity  and 
metal  resistivity)  [16-20,28].  Whereas  the  existing  time-domain 
characterization  methods  are  usually  based  on  the  use  of  the  classical 
RC-  or  RLC-transfer  function  models  [8-15,29].  So,  more  complete 
modelling  method  is  still  necessary  for  the  estimation  of  the  transient 
response  induced  by  the  interconnection  TLs  particularly  for  the  high¬ 
speed  digital-  and/or  mixed-electronic  systems. 

This  statement  motivates  us  to  suggest  and  to  develop  further 
modelling  process  allowing  to  forecast  simply  the  influences  of 
RF / digital  interconnection  circuitries  on  the  signal  quality  via  the 
analysis  of  the  transient  response  behaviours.  To  make  this  paper  more 
understandable,  the  theoretic  analysis  highlighting  the  determination 
of  the  microstrip  TL  high  frequency  parameters  as  the  characteristic 
impedance  and  constant  propagation  from  the  geometrical  (width, 
length,  substrate  height)  and  physical  (dielectric  permittivity  and 
metal  resistivity)  characteristics  is  recalled  in  Section  2.  Then,  the 
calculation  method  enabling  the  extraction  of  the  RLCG-model  per 
unit  length  parameters  is  equally  presented.  Hence,  mathematical 
analysis  explaining  the  determination  of  the  linear  transfer  function 
enabling  to  express  the  interconnection  transient  responses  is  offered 
in  Section  3.  Hence,  discussions  on  validation  results  based  on  the 
comparison  with  standard  electronic  simulation  tool  SPICE  results  are 
the  object  of  Section  4.  Finally,  the  last  Section  is  the  conclusion  of 
the  paper. 
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2.  EXTRACTION  OF  THE  MICROSTRIP  LINE 
RLCG-PARAMETERS 

For  starting  let  us  consider  the  microstrip  TL  with  geometrical  length 
d.  width  w,  metallisation  thickness  t  as  depicted  in  Fig.  2.  One  assumes 
that  the  dielectric  substrate  is  characterized  by  its  height  h.  relative 
permittivity  er  and  loss  constant  tan(<5). 


Figure  2.  Geometrical  representation  of  microstrip  TL. 


The  corresponding  characteristic  impedance  Zc  and  the  propaga¬ 
tion  constant  y(/)  =  a(f)+j  ■  /3(f)  can  be  deduced  from  these  physical 
properties  according  to  the  formulations  reported  in  [16, 17,32]  as  ex¬ 
pressed  in  (1)  and  (3). 


where 


Vo 


In 


h  ■  a[w/h ) 
w 


a(w/h)  =  6  +  (2-7T  —  6)  •  exp 


30.66/A  0,7528 
w  J 


(1) 
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with  T] o  is  the  air  impedance  and  sejj  is  the  effective  permittivity  of 
the  medium.  One  points  out  that  as  reported  in  [32],  by  denoting  c  is 
the  speed  of  light  in  the  vacuum,  the  constant  propagation: 

7  (/)  =  ac(/)  +  ad(f)  +  j^Ver(/),  (3) 

includes  the  metallic  conductor  and  the  dielectric  losses  respectively 
given  by: 
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Figure  3.  TL  RLCG-model. 
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For  the  calculations  including  the  radiating  loss,  more  explicit 
analytical  expressions  of  the  propagation  constant  real  part  or  the  per 
unit  length  loss  constant  a(/)  are  presented  in  [30-33].  Knowing  Zc 
and  7 (/),  one  can  extract  the  equivalent  TL  RLCG-model  with  the 
per  unit-length  parameters  Ru,  Lu.  Cu  and  Gu  as  described  in  Fig.  3. 
As  established  in  [28] ,  the  TL  global  parameters  R  =  Ru  ■  d,  L  =  Lu  ■  d, 
C  =  Cu  ■  d,  and  G  =  Gu  ■  d  can  be  calculated  with  the  following 
expressions: 


R  =  &ky(f)-Zc], 


» [7  (f)-Zc] 


UJ 
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where  9ft  [2]  and  A[z]  are  respectively  the  real  and  imaginary  parts  of  the 
complex  number  z,  and  u>  =  2nf  is  the  operating  angular  frequency. 

By  denoting  p  the  resistivity  of  the  metallic  conductor  constituting 
the  TL  structure,  it  yields  from  the  previous  expressions  that  the  global 
RLCG-parameter  can  be  written  as  follows: 


R  = 


d\/ir  ■  po  ■  p  ■  f 
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(14) 
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G  =  d  ■  tan  (<5)  •  u  ■  C,  (15) 

where  is  the  vacuum  permeability.  By  using  these  parameter 
definitions,  one  proposes  to  establish  the  second  order  transfer  function 
of  the  microstrip  line  before  the  determination  of  the  TL  transient 
responses. 


3.  POLYNOMIAL  MODELLING  OF  THE  MICROSTRIP 
INTERCONNECTION  TL  TRANSFER  FUNCTION 


As  example  of  numerical-digital  interconnection  circuitry  [11,25-27], 
let  us  consider  the  microwave- digital  system  comprised  of  a  TL  driven 
and  loaded  by  logic  gates  represented  in  Fig.  4.  As  can  be  seen  here,  the 
TL  is  characterized  by  its  parameters  ( Zc ,  7)  and  its  physical  length 
d.  One  assumes  that  the  driven  gate  output  behaviour  is  consisted 
of  a  voltage  source  Vs  with  series  impedance  Zs  and  the  load  gate  as 
its  input  impedance  denoted  Zi.  According  to  the  microwave  theory, 
the  transfer  matrix  of  the  interconnection  line  characterized  by  ( Zc ,  7) 
with  geometrical  length  d  is  expressed  as: 


[ttl\ 


cosh  (7  ■  d)  Zc  sinh  (7  ■  d) 
^£0  cosh  (7  -d) 


(16) 


By  considering  the  RLCG-modelling  behaviour  depicted  in  Fig.  3, 
the  overall  interconnect  system  presented  in  Fig.  4  will  be  naturally 
transformed  as  shown  in  Fig.  5. 


gate  source 


interconnection  line 


load  gate 


Figure  4.  Interconnection  TL  driven  and  loaded  by  logic  gates. 

According  to  the  signal  processing  theory,  the  frequency  band  of 
any  digital  signals  is  delimited  by  its  theoretical  bandwidth  where  be¬ 
longs  more  than  85%  of  the  baseband  harmonic  components  [30] .  Thus, 
it  was  established  that  for  the  square  wave  digital  signal  presenting 
rise-/fall-times  tr  the  frequency  band  is  usually  approximately  equal 
to: 

/max  ~  (17) 
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Figure  5.  Interconnection  RLCG-model  driven  by  a  voltage  source 
Vs  with  inner  impedance  Zs  and  loaded  by  Z^. 


It  means  that  to  study  the  analog-digital  systems,  one  can  consider 
its  transfer  function  behaviour  from  DC  to  /max.  More  generally,  one 
can  take  into  account  the  transfer  matrix  of  the  TL  as  its  baseband 
modelling  behaviour  by  exploiting  for  example,  the  equivalent 
polynomial  expression  of  the  matrix  [Ttl]-  This  approximation  can 
be  realized  for  example,  by  applying  the  second  order  MacLaurin 
expanding  of  each  element  of  the  transfer  matrix.  Therefore,  one 
obtains  the  following  expression: 


[Ttl\ 


,  Z(s)-Y{s) 

l-l-  2 

Y(s) 


Z(s) 

i  ,  Z(s)-Y(s)  . 

1  “T  2 


(18) 


with  Z(s)  =  R  +  L  ■  s  and  Y(s)  =  G  +  C  ■  s.  Therefore,  according  to 
the  theory  of  electronic  circuit  and  system,  the  overall  transfer  matrix 
of  the  system  schematized  in  Fig.  5  can  be  assimilated  as: 


[T](s) 


1  Zs 
0  1 
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where 


1  Zs 
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and 
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L  Zjj 


0 
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represents  respectively  the  transfer 


matrices  of  the  series  impedance  Zs  and  the  output  shunt  impedance 
Zl.  From  the  transfer  matrix  [X](s)  expressed  in  (19),  one  can 
demonstrate  easily  that  the  overall  transfer  function  H(s)  can  be 
formulated  as  follows: 


H{s) 


1 

Wi Js) 


1 

ao(s)  +  ai(s)  •  s  +  (s)  ■  s 2  ’ 


(20) 
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where 

a0(s) 

ai(s) 
0,2  (s) 


V  +  R.G  +  l.G.Z.)  .Zl+(1±^G)  Zi  +  ^  {21) 

(R  -  C+L  -  G+2  ■  C  -  Zs)  -  Zl  Ui  (tt-C+L-G)  z  ^ 

2  •( ZL  +  ZS ).  (23) 


To  confirm  the  efficiency  of  this  transfer  function  model,  let  us  see  its 
time  domain  responses  computed  with  Matlab  and  then  compare  with 
the  results  from  the  SPICE-ADS  simulations.  Thus,  we  will  evaluate 
the  examples  of  interconnection  transient  parameters  as  rise- / fall-times 
and  the  50%  propagation  delay. 


4.  VALIDATION  RESULTS 

To  give  evidence  to  the  relevance  of  the  previous  theoretic  concepts, 
let  us  consider  the  arbitrarily  folded  microstrip  interconnection  line 
pictured  in  Fig.  6.  This  microstrip  device  was  simulated  in  SPICE 
schematic  environment  of  standard  microwave  electronic  tool  ADS 
from  Agilent™. 


Figure  6.  3D-design  of  the  simulated  microstrip  interconnection  line. 

This  structure  was  printed  on  the  FR4-substrate  having 
physical  characteristics  and  arbitrary  chosen  geometrical  parameters 
summarized  in  Table  1. 


Table  1.  Physical  characteristics  of  the  tested  microstrip  line. 


w 

h 

t 

£r 

(tan  6) 

500  pm 

625  pm 

35  pm 

4.4 

0.012 
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By  applying  the  RLCG-model  extraction  formulations  introduced 
earlier  in  Section  2,  one  obtains  the  following  per  unit  length  values: 
Ru  =  16.6  0/m,  Lu  =  464.7  nH/rn,  Cu  =  74.7pF/m  and  Gu  = 
5.6mS/m.  With  this  numerical  value,  we  have  implemented  the  above 
transfer  function  into  Matlab  program.  Then,  to  generate  the  transient 
responses  of  the  structure,  a  periodical  trapezoidal  pulse  presenting 
the  temporal  characteristics  (tr:  rise  time,  tf.  fall  time,  Tp:  period, 
tw:  pulse  time  duration)  with  steady  states  Viow  =  OV  and  Vhigh 
normalized  at  IV  addressed  in  Table  2  was  considered  as  the  input. 


Table  2.  Parameters  of  the  considered  input  excitation  signal. 


tr 

TP 

tw 

Vlow 

Vhigh 

100  ps 

lOOps 

1  ns 

400  ps 

OV 

IV 

Since  the  rise- /fall-time  of  the  considered  input  signal  is  about 
lOOps,  according  to  (10),  the  maximum  frequency  required  for  the 
proposed  polynomial  model  response  have  been  estimated  of  about 
3.5  GHz.  We  verify  that  it  enables  to  generate  transient  response  taking 
into  account  the  propagating  TEM-modes  through  the  under  test 
structure  which  is  compliant  with  the  second  order  transfer  function 
expressed  in  (20).  In  order  to  demonstrate  the  adaptability  of  the 
proposed  modelling  method  for  various  types  of  the  load  impedance,  let 
us  analyse  the  computed  results  first,  for  the  resistive  load  (Zl  =  Rl) 
and  then,  for  the  capacitive  load  (Zl  =  l/( Cl  ■  p))  in  the  next 
paragraphs.  As  aforementioned,  the  comparative  results  are  referred 
with  SPICE  simulations  run  in  ADS  transient  simulations  regarding 
the  physical  parameters  summed  up  in  above  Table  1. 

4.1.  Transient  Responses  for  Zl  =  Rl 

It  is  interesting  to  note  that  the  static  gain  of  the  overall  understudy 
system  (see  Fig.  4)  is  equal  to  H( 0)  =  Rl/Rs  when  the  source  and  load 
impedances  are  purely  resistive.  In  the  present  case,  we  have  changed 
the  value  of  the  resistive  load  Rl  =  {10  fi,  500  H}  and  the  length  of  the 
microstrip  line  d  =  {8.5  mm,  13.5  mm}.  So,  one  gets  the  comparative 
results  plotted  in  Fig.  7.  The  proposed  model  responses  are  plotted  in 
full  lines  and  the  ADS  simulation  results  are  in  dashed  lines. 

One  can  see  that  the  transient  response  generated  from  the 
proposed  model  agrees  very  well  with  the  SPICE  simulations.  One 
observes  that  through  the  microstrip  interconnection  the  pulse  signal 
is  significantly  degraded.  Here,  with  the  two  computation  processes, 
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one  evaluates  rise-/fall-times  of  about  tr  ~  0.14  ns  for  Rl  =  50014  and 
tr  ~  0.43  ns  for  Rl  =  1014.  Then,  one  evaluates  rise-/fall-time  relative 
errors  of  about  1%. 

4.2.  Transient  Responses  for  Zl  =  1  /{Cl  •  p ) 

In  this  case,  Zl  is  assumed  as  a  capacitive  load  C^.  In  order  to  show 
the  interaction  between  these  load  values  and  the  interconnect  line, 
we  have  changed  the  values  of  the  capacitive  load  Cl  =  {lpF,  10 pF, 
30  pF}  and  the  lengths  of  the  microstrip  line  d  =  {8.5  mm,  13.5  mm} 
to  demonstrate  the  flexibility  of  the  proposed  modelling  process.  As 
consequence,  one  gets  the  comparative  results  displayed  in  Fig.  8  below. 
The  proposed  model  responses  are  plotted  in  full  lines  and  the  ADS 
simulation  results  are  in  dashed  lines. 


time  (ns) 

(b) 

Figure  7.  Transient  responses  of  the  microstrip  structure  shown 
in  Fig.  6  for  different  values  of  the  load  gate  resistance  Rl  for  (a) 
d  =  8.5  mm  and  (b)  d  =  13.5  mm. 


Progress  In  Electromagnetics  Research,  Vol.  112,  2011 


193 


time  (ns) 


(a) 


(b) 


Figure  8.  Transient  responses  of  the  microstrip  structure  shown  in 
Fig.  6  for  different  values  of  the  capacitive  load  Cl-  (a)  d  =  8.5mm 
and  (b)  d  =  13.5  mm. 

One  can  see  that  the  output  is  completely  distorted  for  the 
capacitive  load  C l  higher  than  1  pF.  One  evaluates  here  a  50% 
propagation  delay  of  about  T„50%  =  {80  ps,  240  ps,  440  ps}  respectively 
for  Cl  =  {lpF,  10 pF,  30  pF}. 

For  the  both  examples  of  the  performed  validation  results,  one 
points  out  that  the  computation  time  with  the  proposed  method  is 
in  order  of  tens  microseconds.  The  presented  method  is  practically 
simpler  than  most  of  3-D  EM  solver  ones  as  the  example  of  the  TD 
integral  equation  presented  in  [34]. 

5.  CONCLUSION 

A  transfer  function  modelling  method  of  the  millimetre  interconnection 
for  the  high-speed  integrity  analysis  in  the  RF-/digital  PCB  is 
investigated.  The  established  model  is  based  on  the  exploitation  of 
the  electrical  RLCG-model  transfer  matrix  which  was  assumed  as  its 
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second  order  polynomial  linear  model. 

To  evidence  the  functionality  of  the  introduced  method,  a 
microstrip  interconnection  TL  driven  and  loaded  by  the  logic  gates 
was  considered  and  analysed.  It  was  described  that  the  per  unit 
length  equivalent  model  parameters  were  determined  from  the  TL 
physical  and  the  geometrical  parameters  as  the  substrate  permittivity 
and  height,  and  the  line  width  and  length.  Then,  the  equivalent 
transfer  function  linear  model  of  the  overall  system  composed  of  a  TL 
combined  with  the  logic  gates  is  mathematically  established.  Thus, 
for  the  testing  process,  it  was  excited  with  a  periodical  trapezoidal 
pulse  voltage  with  2  Gbits/s  rate.  So,  good  agreement  between  the 
transient  responses  from  the  established  transfer  function  calculated 
with  Matlab  and  those  from  SPICE  simulations  of  the  overall  structure 
was  found.  It  was  shown  that  the  microstrip  interconnection  responses 
are  completely  degraded.  Then,  after  the  sweep  of  the  load  impedance 
values  assigned  as  resistance  and  capacitance  loads,  one  observes  also 
that  the  calculated  transient  results  remain  well-correlated  to  the 
simulations. 

One  emphasizes  that  the  design  and  the  simulation  of  this  type 
of  microwave  digital  interconnections  become  very  difficult  to  carry 
out  when  the  circuit  is  composed  of  thousands  of  logic  gates.  For  this 
reason,  it  seems  important  to  exploit  the  proposed  modelling  method 
for  the  prediction  the  analogue-digital  signal  behaviours  along  the 
interconnections.  Thanks  to  the  accuracy  and  the  computation  time 
gain,  we  think  that  the  proposed  method  can  be  a  good  candidate  for 
the  modelling  of  the  complex  structure  of  the  interconnection  circuitry 
in  the  high  density  integrated  circuit  as  the  clock  distribution  networks. 

In  the  continuation  of  this  work,  one  plans  to  improve  this  method 
for  the  estimation  of  the  RF-microwave  and  digital  electronic  device 
interconnection  effects  by  taking  into  account  the  non-uniformity  of 
the  TL  as  introduced  in  [35]  and  the  eventual  crosstalk  with  the 
neighbourhood  TLs  [36].  Then,  we  would  employ  the  proposed  method 
to  the  development  of  the  equalization  technique  with  the  NGD  circuit 
for  the  reduction  of  the  correction  of  the  digital  signal  degradation  with 
MMIC  and  digital  integrated  systems. 
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